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Abstract 

Time reversal invariance can be summarised as follows: no difference can be measured if a 
sequence of events is run forward or backward in time. Because price time series are domi- 
nated by a randomness that hides possible structures and orders, the existence of time reversal 
invariance requires care to be investigated. Different statistics are constructed with the property 
to be zero for time series which are time reversal invariant; they all show that high-frequency 
empirical foreign exchange prices are not invariant. The same statistics are applied to math- 
ematical processes that should mimic empirical prices. Monte Carlo simulations show that 
only some ARCH processes with a multi-timescales structure can reproduce the empirical find- 
ings. A GARCH(1,1) process can only reproduce some asymmetry. On the other hand, all the 
stochastic volatility type processes are time reversal invariant. This clear difference related to 
the process structures gives some strong selection criterion for processes. 



Keyword: Time reversal symmetry, ARCH processes, stochastic volatility processes. 
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1 Introduction 



Time reversal invariance (TRI) is a very important concept in science. The idea can be sum- 
marised as follows: when a sequence of events is viewed starting from the end, namely with 
the arrow of time reversed, is it possible to measure a difference compared to the normal time 
ordering? A rigorous formulation of time reversal invariance is that the transformation t —> —t 
is an exact symmetry of the system under consideration. The basic laws of physics are time 
reversal invariant (Newton equation for mechanic, Maxwell equations for electromagnetism, 
Einstein equation for general relativity, Dirac equation for quantum mechanics, etc.), but the 
macroscopic world is clearly not time reversal invariant. This paradox was solved by thermo- 
dynamics and the increase of entropy. The same question can be asked about finance, namely 
if a time series of prices originating in the financial market is time reversed, can we "see" the 
difference? Contrarily to a Buster Keaton movie, at the level of prices or returns, it is indeed 
very difficult to notice a difference because financial time series are dominated by a randomness 
that hides possible structures and orders. Therefore, the appropriate formulation of the question 
is whether statistics can show the presence, or absence, of time reversal invariance. 



[Ramsey andRothman, 1988, Ramsey and Rothman, 1996 1 have already addressed this ques 



tion in an economic framework. Their idea is to search for differences between up and down 
moves over long time horizons, using yearly economic indicators. This behaviour is typically 
related to business cycles, for example with long slow rises followed by abrupt decreases. The 
estimator that is used is given by E [ r 2 (t) r(t + kdt) ] where bt is the time increment of the 
time series, r the return and k ^ an integer index. As the amount of long term economic 
data is fairly small, such studies need to rely on a carefull analysis of the statistical prop- 



erties of the indicators | Ramsey andRothman, 1996 1. Other indicators have been proposed 



[ |Chen et al., 2000| , Fong, 2003 1. The salient results is that most of the economic time series are 



not time reversal invariant, but for a small fraction of them, the null hypothesis of time reversal 
invariance cannot be rejected. 

This paper takes a different angle to study TRI in empirical data by using high-frequency for- 
eign exchange time series. The focus is to study various statistics related to the volatility, for 
time horizons ranging from 3 minutes to 3 months. For foreign exchange rates, a symmetry 
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between the exchanged and expressed currencies is plausible (at least for major free floating 
currencies). This symmetry occurs because a FX rate is a conversion factor between two nu- 
meraires, and not the price of a security expressed in a numeraire (like for an equity price or 
of a bond price). Under the exchange of the currencies, an exchange rate is transformed by 
price — > 1 /price, and the logarithmic returns by r — > — r. Notice that reversing the time induces 
the same transformation r — > —r and the reverse ordering of the time series. Therefore, if the 
exchange of currencies is an exact symmetry, then all statistics that are odc£ in the returns are 
zero, like for example E [ r 2 (t) r(t + kdt) ] . The same argument implies that the return proba- 
bility distribution in even, that is, p(r) = p(—r). We have checked below the empirical validity 
of this argument for pir), and found that it is likely incorrect. Even if ultimately incorrect, the 
argument points to a small asymmetry for such term. In order to have a better signature of time 
irreversibility, a better track for foreign exchange data is to search for estimators that are even 
in the returns, but sensitive to time reversal. 



A general definition of TRI is given for example in JChen et al., 2000[ |. In particular, they show 
that the distribution of the returns must be symmetric for a TRI series, and they construct a test 
based on this property using the characteristic function of the pdf. A more general construction 
along this line is as follow. Take a quantity a such that a — > a under the time reversal trans- 
formation. For example a can be the logarithmic price; in this paper a is a volatility estimator 
computed from the returns. Then, the quantity Ao(/) = o(t + dt) — a(t) is odd under the time 
reversal transformation: Ao — > — Ao. If the series is TRI, then the distribution of Aa must be 
even. Another idea for testing TRI is based on covariance or correlation between two quanti- 
ties x and y. Essentially, the quantity E [ x(t) y(t + dt) ] should be equal to E[x(t + 8t) y(t) ] if 
the series is TRI, and a test can be constructed on the difference between these two quantities. 
Because of the symmetry between the arguments of the covariance or the correlation, the two 



quantities must be different in order to have a non trivial test. [ [Ramsey and Rothman, 1988) use 
r and r 2 ; this paper uses volatilities a with different parameters. 

This paper presents three different statistical estimators sensitive to TRI. These statistics are 
based on various estimator of the volatility (i.e. even in the returns), and are essentially measur- 
ing time reversal for volatilities computed with information before t and after t. They all show 



'An even function is such that f(—r) = /('*), an odd function such that f(—r) = —f(r) 
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that empirical data are clearly not time reversal invariant. 

Using the same statistical tool, various processes that should mimic empirical data are inves- 
tigated. For example, a simple Gaussian random walk is time reversal invariant, and the three 
statistics are zero. Much more interesting are the processes from the ARCH family and from 
the stochastic volatility family. Some of the ARCH type processes can reproduce the empirical 
figures, but all the simulated stochastic volatility processes are time reversal invariant. These 
results demonstrate that the last family of processes cannot describe some stylized facts of fi- 
nancial time series. It is an important results as TRI statistics allow us to select between models 
that are structurally very different. The usual method for model selection is to nest the pro- 
cesses and to show that the corresponding parameter is significantly different from zero. But 
this method is not generally possible for processes that are so widely different as ARCH and 
stochastic volatility. Besides, the systematic comparison of the results between empirical data 
and processes allow us to glimpse the origin of time asymmetry. This discussion is presented 
just before the conclusion. 



2 Data and notations 

The empirical data used in the empirical study originate from the foreign exchange market. 
The high-frequency tick-by-tick quotes are used to compute a continuous "business time scale" 



[Breymann et al., 2000 1 in order to deseasonalize the strong daily and weekly patterns present 
in such data. The prices are sampled each 3 minutes in business time in order to obtain de- 
seasonalized homogeneous time series. An efficient deaseasonalisation procedure is crucial for 
the presented computations, as the use of a 3 minutes regular sampling in physical time (after 
the week-ends have been removed) would mainly show the strong daily pattern related to the 
opening and closing of the various market around the world. Such strong seasonality would 
hide other interesting stylized facts. The author is gratefull to Olsen & Associates in Zurich for 
providing the deseasonalized data. The in-sample used in this study starts January 1, 1990 and 
ends July 1, 2001. The year 1989 is used to build-up the computations. 

Our notations are as follows: r denotes a time series and r(t) the value of the time series at time 
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t. The parameters are denoted between bracket, like for the time series r[dt r ] or for the value 
r[dt r ] (t) at time t. The subscript are used to denote different estimators, for example (J/, and a r . 
The raw data are given by the logarithmic price time series x(t) defined between a start time 
t s and an end time t e . The time reversal transformation corresponds to changing x(t s + AT) by 
x(t e — AT), which is denoted informaly as x(t) — > x(—t).. The return r over a time interval ht r 
is defined by r[5t r ] (t) = x(t) - x(t - dt r ) . 

The historical volatility Oh is 

a 2 h [8t a M(t)= 1 -^- £ r 2 [St,](t'). 

° h n t-8t G +bt r <t'<t 

Essentially, 0/ 7 (7) measures the fluctuation of the prices at the scale 8t r , in the time interval [t — 
8t a ,t]. This definition includes information only in the past of t. The ratio lyear/67 r annualises 
the volatility, and n is the number of terms in the sum over t' . In the actual evaluation, the sum 
over t' is carried over each point on the 3 minutes time grid. 

The realized volatility a,- is essentially the same definition, but using information only in the 
future of t 

O r [8f , St r ] (t) = £ ? [Str] (0 = O h [& , Str] (t + &b) • 

°' r n t+8t r <t'<t+8t a 

Notice that a definition of volatility depends on two time intervals, namely the time interval 
dt r over which the returns are computed (also called the granularity) and the time interval bt a 
over which the returns variance is computed. In order to have a good volatility estimator, the 
ratio dt a /§t r should be large enough, say dt a /& r > 10. Except for the third time asymmetry 
estimator, this ratio is fixed to dt a /dt r = 24. 

3 Empirical time reversal statistics 

The changes of the volatility are measured by the volatility increment 

Ao = G, -G/,. 
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With a time reversal transformation, the volatility increment changes by Ao(t) — > — Ao(t). 
The probability density p(Ao) of the volatility increment can be estimated and its asymme- 
try a p (Ao) = p(Ao) — p(—Ao) gives a measure of the time irreversibility. This asymmetry is 
a quantitative measure of the following intuitive perception of the price dynamics. A shock on 
the market (for example due to the arrival of an important piece of news) produces a sudden 
increase of volatility, followed by a slow relaxation toward a normal level of volatility. For 
the time reversed series, this corresponds to a slow volatility increase followed by a sudden 
return to the normal and the distribution of volatility increments is p(— Aa). The asymmetry 
a p (Ao) measures the asymmetry in the dynamics between the original and time reversed series. 
Figured] shows the probability distribution p(Ao). The probability is estimated by binning the 
empirical data using a linear interpolation in a non uniform grid. The points on the sampling 
grid are chosen so that a roughly equal number of values fall into each bins. 

The pdf appears to be symmetric at first glance. It shows that our intuitive perception as de- 
scribed above is exaggerated. Yet, a detailed examination of p(Ac) reveals the expected asym- 
metry around Aa ~ 0.02. Figure [2] displays the asymmetry a p (Ao), and a fairly consistent 
symmetry breaking pattern for various empirical time series is observed. The negative values 
for Ao < 0.05 corresponds to the "return to the normal" or to a larger probability for small 
negative volatility increments. There is a corresponding larger probability for large volatility 
increments Aa (that is the arrival of news or shocks), that translate into positive values for a p . 
Because the volatility is stationary, a simple empirical first moment like (Aa) converges toward 
zero with the inverse of the sample size. This shows that more complex statistics should be used 
to reveal this asymmetry. 

The second statistics involves the correlation of volatilities at various time horizons. The fol- 
lowing correlation is investigated in [ |Zumbach and Lynch, 200 1[ Lynch and Zumbach, 2003 1: 



p a (dt ai dt' a ) = p (a H [8t a ,& a /2A](t),a r [K> K/ 24 } (0) 

where on the right hand side p(x,y) is the usual linear correlation between two time series x 
and y. Essentially, this quantity measures the dependency between past and future volatilities, 
at the respective time horizons dt a and 8t' a . This correlation proves to be a very powerful tool 
to investigate the structure of the underlying market; in particular, it shows that the markets are 
heterogeneous as structures can be observed at the natural human time horizons (intra-day, day, 
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Figure 1: Probability density for the volatility increment Ao, for the foreign exchanges 
CHF/USD, DKK/USD, JPY/USD, USD/GBP and for gold XAU/USD. The time horizons are 
dt a =1 day, ht r = dt a /24. Similar figures are obtained for other time horizons. 
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Figure 2: Asymmetry of the probability density for the volatility increment Aa. The parameters 
and time series are as for fig. [7J 

week and month). Overall, this correlation is asymmetric with respect to the exchange of dt a 
and bt' a . As this exchange is directly related to the time reversal symmetry, a historical versus 
realized volatility correlation asymmetry is defined by 

a a (dt a , &o) = p (& , &„) - p (54 Sf ). 

and a a ~ indicates time reversal invariance. The computation of a a for the above empirical 
time series reveals a fairly consistent pattern, with a maximum of order 6 to 12% for dt a ~ 1 
week and bt' a ~ 6 hours. The natural representation for a a (dt a ,dt' a ) is in a two dimensional 
plane, but the main behaviour can be capture in a 1 dimensional cut, as given in fig. [3] The 
empirical data show a distinct and fairly consistent pattern indicating a clear asymmetry with 
respect to time reversal. 

The third statistics involves correlations of past and future volatilities at the same time hori- 
zon 8f , but different granularities 8t r and ht' r respectively. For a given ht Q , the "granularity" 
dependency is defined by 

Pgr(5f r ,5r;) = p(a h [8t c ,8t r }(t),a r [dt a ,8t' r ](t)) 
7 




02:24 03:36 08:24 15:36 1d08:24 2d 16:48 5d 09:36 10.75d 21.55d 

volatility time horizon 

Figure 3: The measure of asymmetry a (8? ,8^). The parameter dt' a is given on the horizontal 
axis, the value of dt a is the symmetric through the vertical axis (for example for dt' a = 8h24 
corresponds dt G = 5 day 9h36). 
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00:03 00:06 00:12 00:24 00:48 01:36 03:12 06:24 12:48 

volatility time horizon 

Figure 4: The measure of asymmetry a gr (8t r , 8t' r ). The parameters are: dt a = 2 9 • 3 minutes = 
1 days lh36; dt' r = 2" • 3 minutes is given on the horizontal axis; dt r = 2 8 ~" ■ 3 minutes is the 
symmetric value through the middle point (for example to ht' r = 24 minutes corresponds dt r = 
lh36). 
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with bt r < dt a , §t' r < dt a . In [ |Dacorogna et al., 1998[|Dacorogna et al., 2001) a similar statistics 
was introduced, showing the asymmetry between fine and coarse volatilities. As for the second 
statistics, the exchange of dt r and dt' r is related to the time reversal symmetry, and a similar 
measure of granularity asymmetry is defined by a gI (dt n dt' r ) = p gT (8t ri 8t' r ) — p gr (8^,&,-). The 
computation of a gI for empirical time series at dt a ~ 2 days shows a consistent and systematic 
asymmetry, with values in the range of 12 to 18%. The natural representation is in the two 
dimensional space (&t r , 8t' r ), and a one dimensional cut is displayed on figured The asymmetry 
pattern is clear and consistent between the empirical data. Therefore, the three measures of 
asymmetry deliver a consistent message, namely foreign exchange time series are not time 
reversal invariant, and the asymmetry is quantatively small. 



4 TRI in theoretical processes 

A similar study can be conducted for theoretical processes, using Monte Carlo simulations. 
The simplest process is a Gaussian random walk, that is exactly time reversal invariant (the 
proof follows from the independence of the increments that allows to reorder terms under the 
expectations). The more interesting processes include heteroscedasticity, either with an ARCH 
form or with a stochastic volatility term (see e.g. QPoon, 2005] ] for a recent general reference). 



In an ARCH process, volatility is a function of previous returns. This function can be fairly 
general and include for example multiple time horizons, or various powers of the returns (r 2 , 
\r\, ln( | r | ) , ...). An investigation among several ARCH processes shows that they are not 
time reversal invariant. The simplest GARCH(1,1) process [ [Engle, 1982 |Bollerslev, 1986 



Engle and Bollerslev, 1986 1 exhibits an asymmetry according to the measures 1 and 2, but the 
last measure displays no asymmetry (see figs. [2] [3] and H]). The multiple time horizon processes 
introduced in [|Zumbach and Lynch, 200 1[ |Lynch and Zumbach, 2003[ |Zumbach, 2004] can re- 



produce quantitatively the above three measures of time irreversibility. A key ingredient to put 
in these multiscale processes is that the return time horizons dt r must increase with the volatility 
time horizon dt a . This corresponds to the intuition that short term intra-day traders use tick-by- 
tick data, whereas long term fund managers use daily data. When dt r is kept at the process time 
increment, the time reversal asymmetry is too small, or even zero for the third measure. Taking 
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a volatility granularity of order 5t a = 24 ■ dt r gives roughly the correct quantitative time reversal 
asymmetry in these multiple components processes. On the other extreme, taking dt a = dt r 
produces too much asymmetry. Clearly, the process parameters can be chosen in order to reach 
a better numerical agreement with the empirical values (no such optimisation has been done for 
this work). 

In a stochastic volatility (SV) process, the volatility is an independent process with its own 



source of randomness (see e.g. [Shephard, 2003 1 for a general reference). The return is a 



"slave" process dependent on the volatility, and there is no feed-back from the return to the 
volatility. The volatility process can be fairly simple, as in an exponential SV process or as 
in the Heston process. In the exponential SV process, the logarithm of the volatility follows 
an Ornstein-Uhlenbeck process. In the Heston process, the volatility follow a random walk 
with mean revertion, but the stochastic term is multiplied by y/a so that a stays positive. Both 
processes can also be extended to include multiple time horizons in order to induce a volatility 
cascade over multiple time scales. Two such models with a long memory volatility cascade are 
label in the graphs and tables as "LM stoch.vol." and "LM Heston". Yet, all these SV processes 
show time reversal invariance (up to the statistical errors of the numerical simulation, see figs.[2l 
[3] and©. 

A family of models with a similar structure is the regime switching processes. In this case, 
the volatility is given by another independent process with an integer value that gives the "state 
of the world" (like quiet, excited, shocked, etc.). To each index value corresponds a volatil- 
ity, and the return follows a simple random walk with the corresponding volatility. The state 
of the world is an independent process, without feed-back from the return, and follows a sim- 
ple Markov process with probability Pr(z — > j) to jump from state i to state j. The transition 
probabilities Pr(z — > j) are strongly constrained to give realistic distribution for the volatility. 
With asymmetric transition probabilities Pr(z — > j) ^ Pr(j — > z), the process is not time reversal 
invariant according to the measures a p and a a , as is visible on fig. |5j Yet, the asymmetry is 
much smaller than the empirical observed values, and it is not possible to modify the transition 
probabilities to get simultaneously large asymmetry and realistic pdf. Finally, the coarse/fine 
graining measure of asymmetry a gr for the regime switching process is compatible with zero, 
up to Monte Carlo statistical fluctuations. Because of the lack of feed-back of the return on the 
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Figure 5: Asymmetry of the probability density for the volatility increment Ao for a few pro- 
cesses. 

volatility, one can understand intuitively that the stochastic volatility and regime switching can- 
not include an asymmetry in the third measure of invariance a gr using the volatility granularity. 



5 Test statististics 



So far, the focus was to construct statistics sensitive to time reversal invariance, to understand 
what they are measuring, and to show the differences between empirical time series and various 
processes. Yet, definitive conclusions require to use rigorous test statistics. Essentially three 
methods can be followed in order to obtain distribution information on a statistics: analytical, 
bootstrap and Monte Carlo simulations. In order to select the most appropriate method for the 
asymmetry statistics, one should keep in mind that our statistics are based on volatilities and 
that volatilities have slowly decaying correlations (like a power law). The bootstrap method 
is appropriate only when the data are independent, or possibly when the dependencies decay 
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exponentially fast. Clearly, the slow decay of the correlations rules out a bootstrap approach. 
Similarly, the analytical approaches are using a convergence toward a limit law, like a Gaussian. 
The convergence needs to be sufficiently fast, and similar conditions are imposed on the data as 
for the bootstrap method. Moreover, the asymmetry statistics are not based on simple moments, 
but on probability distributions and correlations. Both problems makes the analytical approach 
fairly difficult. Therefore, we have used Monte Carlo simulations to compute numerically the 
distributions for the asymmetry statistics and the related p-values. The computed p-values are 
the probability that the statistics is lower or equal to zero. For a symmetic distribution, we 
expect a p-value around 0.5. 

In order to obtain the distribution of the asymmetry statistics for a given process, 200 Monte 
Carlo simulations are performed. For each simulation, the process is simulated with a time 
increment of 3 minutes and for a time length identical to the available empirical data, namely 
11.5 years. At the end of the simulation, asymmetry statistics are computed. This is repeated 
200 times, and the empirical distributions for the asymmetry statistics are computed. The cu- 
mulative probability that a statistic is smaller or equal to zero can be easily estimated, and gives 
the desired p value. This approach is simple and gives finite sample information, but relies on 
the process that should reproduce sufficiently well the empirical data. 

The distribution statistics can be obtained for given values of the arguments, say for example 
a p (Ac) for a given Ag. Yet, the sensitivity of the TRI test can be enhanced by simple integration 
of the above statistics (similarly to a portmanteau statistics for the correlation). The volatility 
increment asymmetry is defined by 



where Ao^ are the points on the histogram sampling grid and where n is the number of terms in 
the sum. In the computations below, the upper bound is Aa = 0.06. For the historical-realized 
volatility correlation asymmetry, two possibilities are 



Ap(Ac) 



1 



(1) 



n 



0<Aa k <Aa 




(2) 



(3) 
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mean 


stdDev 


p-value 


CHF/USD 


0.64 






DKK/USD 


0.39 






JPY/USD 


0.65 






USD/GBP 


0.95 






XAU/USD 


0.63 






GARCH(U) 


0.19 


0.083 


0.014 


LM-Aff-Agg-ARCH 


0.26 


0.082 


0.015 


Mkt-Aff-Agg-ARCH 


0.39 


0.1 


0.0 


exp stoch.vol. 


0.062 


0.12 


0.35 


exp LM stoch.vol. 


-0.02 


0.09 


0.48 


Heston 


0.008 


0.075 


0.36 


LM Heston 


0.02 


0.1 


0.50 


Regime Switching 


0.23 


0.11 


0.015 



Table 1: The total measure A p for the asymmetry of the probability density of Ao at a time 
horizon of 8t a = 1 day. The columns give respectively the mean, standard deviation and p- 
value for A p . 

The first quantity computes the total asymmetry in the historical-realized volatility, the second 
one the asymmetry along the cut used for figure [3j In practice, the two statistics give very 
similar results. Table [2] reports the p-values for A a cut as the values are related to figure [3] For 
the volatility graining asymmetry, the same two integrated measures can be constructed. Again 
to stay close to figure H the statistics for 

/■8/ max 

A gr ,cut= / dt r ag r(St r ,f(8t r )) (4) 

J 8'min 

are reported. 

The main statistics are reported in the table \T\ to [3j For the historical versus realized volatil- 
ity correlation asymmetry, the probability distributions for the processes are given in figure [6J 
The p-statistics confirm fully the above picture: the empirical data are compatible only with 
multicomponent ARCH processes. Notice that the process parameters have not been optimized 
to reproduce quantitatively the empirical asymmetry. A better agreement can be obtained for 
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mean 


stdDev 


p- value 


CHF/USD 


0.1 






DKK/USD 


0.09 






JPY/USD 


0.05 






USD/GBP 


0.09 






XAU/USD 


0.08 






GARCH(1,1) 


0.007 


0.02 


0.42 


LM-Aff-Agg-ARCH 


0.14 


0.01 


0.0 


Mkt-Aff-Agg-ARCH 


0.15 


0.01 


0.0 


exp stoch.vol. 


0.0009 


0.01 


0.50 


exp LM stoch.vol. 


-0.0016 


0.01 


0.52 


Heston 


0.0005 


0.007 


0.49 


LM Heston 


-0.001 


0.01 


0.54 


Regime Switching 


0.004 


0.01 


0.34 



Table 2: The integrated measure of asymmetry for the historical/realized volatility correlations 
^o,cut along a cut in the (oY a , $t' a ) plane. 
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Figure 6: The probability density for the asymmetry measure A a ^ cut . 
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mean 


stdDev 


n-valne 


CHF/USD 


0.13 






DKK/USD 


0.1 






JPY/USD 


0.095 






USD/GBP 


0.1 






XAU/USD 


0.08 






GARCH(U) 


0.0008 


0.007 


0.48 


LM-Aff-Agg-ARCH 


0.1 


0.01 


0.0 


Mkt-Aff-Agg-ARCH 


0.11 


0.01 


0.0 


exp stoch.vol. 


-0.00007 


0.008 


0.58 


exp LM stoch.vol. 


-0.0007 


0.006 


0.54 


Heston 


0.00007 


0.01 


0.52 


LM Heston 


0.0003 


0.007 


0.59 


Regime Switching 


0.0002 


0.01 


0.46 



Table 3: The integrated measure of asymmetry for the volatility graining correlation A gr , cut 
along a cut in the (8? r , 8t' r ) plane. 
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mean 


stdDev 


p- value 


CHF/USD 


-0.21 






DKK/USD 


-0.23 






JPY/USD 


-0.11 






USD/GBP 


-0.08 






XAU/USD 


0.21 






GARCH(1,1) 


-0.009 


0.07 


0.54 


LM-Aff-Agg-ARCH 


-0.005 


0.07 


0.45 


Mkt-Aff-Agg-ARCH 


-0.0009 


0.09 


0.44 


exp stoch.vol. 


-0.002 


0.08 


0.42 


exp LM stoch.vol. 


0.0015 


0.08 


0.54 


Heston 


-0.002 


0.075 


0.50 


LM Heston 


0.002 


0.08 


0.55 


Regime Switching 


-0.005 


0.08 


0.59 



Table 4: The total measure of the asymmetry of the probability density of r at a time horizon 
of 1 day. The columns give respectively the mean, standard deviation and p-value. 

the multicomponent ARCH processes as they contain parameters that control the asymmetry. 
This is not possible for other stochastic volatility processes, say for example with respect to 
the volatility graining asymmetry: regardless of the parameter values, the p-values stay close to 
1/2. Notice also that a simple GARCH(1,1) process can reproduce the asymmetry for a p , but 
the asymmetry for a a is too small and the asymmetry for a gr is not reproduced at all. 

6 Possible origins of time irreversibility 

With the help of both the empirical analysis and the process structures, we can speculate on 
the origin of the time asymmetry. The systematic study on the processes indicates that the 
time direction is set by the feed-back loop of the price changes on themselves through the 
intermediary of the price volatility, as captured by an ARCH process. Moreover, to obtain the 
measures of asymmetry with the right magnitudes, multiple time scales must be used in an 
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ARCH structure. This indicates that the interplay between the different time horizons as well 
as the feed back loop through the price are enough to create the observed time asymmetry. This 
argument emphasises the role of the prices as the only vector of information between market 
participants. By contrast, other "hidden" variables such as an independent stochastic volatility 
or other market states are irrelevant with respect to TRI. The picture that emerges is of markets 
segmented along the time horizons of the market participants. A more detailed analysis of the 
historical versus realized volatilities points in the same direction, with a cascade from short time 



horizons to long horizons [Zumbach and Lynch, 200 1[ Lynch and Zumbach, 2003 1 



Another relevant topic with respect to TRI is the connection between news and market be- 
haviour. Clearly, financial markets are open systems, driven by various external sources of 
information, like exchange of goods, politics, central banks, etc... Even if intuitively clear, the 
connection between "news" and financial markets is very difficult to establish from a purely 
empirical point. The two fundamental problems are the economic quantification of a "string 
of text", and the discounting by market participants of the "expected" portion of a news. With 
respect to TRI, the important question is whether the irreversibly is of external origin (i.e. from 
the news) or endogenous (i.e. created by the market participants or the trading rules). As such, 
this question is very difficult to investigate empirically, but the present study offers some indi- 
rect evidence. In the equations for a process, there is at least one source of randomness, for 
example in the return equation r(t) = a(t) e(t). Possibly, a process can include several sources 
of randomness, like in a stochastic volatility process. Essentially, the random variables e(t) cap- 
ture two different phenomenon: the trades of the market participants (for example through an 
order queue), and the influence of the external world (i.e. the news). The canonical hypothesis 
for the random processes £ is to assume an iid distribution (and independence at a given time 
with the other process variables like price, return, volatility, etc ...). A debatable assumption 
is that of time series independence: it is plausible that the news are serially correlated, as an 
important piece of information is likely to be followed by more information. But this is dif- 
ficult to establish directly, again because news are given by texts and not by numbers. In the 
investigation of processes, we follow the canonical hypothesis of source of randomness that are 
iid. Within this hypothesis, two processes are particularly interesting with respect to TRI. In a 
regime switching process, the "state of the world" i(t) exhibits irreversibility, as controlled by 
the transition probabilities Pr(z — > j), and i(t) has non trivial serial correlations. For this pro- 



18 



cess, a parallel can be established with the serial correlation of the news. Yet, it is not possible 
with a regime switching process to reproduce the observed time asymmetry. On the other hand, 
within the framework of iid source of randomness, the multiscale GARCH processes can pro- 
vide a quantitatively correct description of the time irreversibility observed in the empirical time 
series. Together, both arguments indicate that the time irreversibility is mostly of endogenous 
origin, namely created by the interactions between market participants, and not from external 
origin. 

In a physical system, a common origin of time irreversibility is in friction or dissipation, and 
in the increase of entropy. A direct analogy of these quantities in a financial system should be 
taken with care. For example, the equivalent of friction can be thought as the bid-ask spread, 
namely the cost for buying, followed by selling directly afterward. However, in a liquid markets, 
and particularly in modern electronic trading exchanges, the market orders and limit orders 
play a very symmetric role. Using an argument based on the optimal choice of the market 



agents between both kind of orders, the authors in pWyart et al., 2006[ show that the profit of a 
systematic market making trading strategy should be close to zero. Indirectly, this shows that 
the raw cost of trading represents a very small dissipation in the system. In the same direction, 
all the simulated processes do not include the cost of transactions, yet all the ARCH processes 
are not TRI. Both arguments show that the observed asymmetry in time does not originate in 
the bid- ask spread. 

Notice that the above discussion on the origin of time irreversibility presents only indirect ar- 
guments, as only an investigation of the market microstructure and of the market participant 
decision processes can give direct evidence. Yet, such studies are clearly very difficult, and 
even more on large decentralised markets as foreign exchanges. 

All the measures of asymmetry used so far are even in the returns, and essentially related to 
volatilities. The asymmetry in the return distribution is another measure of time asymmetry, 
which is odd in the returns. As discussed in the introduction, this symmetry is also related 
to the exchange of currencies for the foreign exchange time series. The empirical results for 
the return distribution asymmetry are reported in table |4] All the processes have rigorously 
symmetric return distributions (because the returns are proportional to the residuals, which 
have a symmetric distributions). The results for the processes are in clear agreement with this 
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symmetry, with p-values close to 1/2. Interestingly, the asymmetry for the empirical data is 
much larger than for any processes (roughly at 1 a to 2.5 a from 0). This indicates that the 
argument on the exchange of currencies in a FX rate, as given in the introduction, is likely not 
correct. Clearly, more extensive investigations are needed to clarify the empirical facts and to 
develop processes that can reproduce the observed return asymmetry in the data. 



7 Conclusions 



To conclude, the systematic study of the time reversal invariance in finance proves to be a very 
powerful investigation tool. The empirical time series are clearly not time reversal invariant, 
according to three different measures. Although this fact is not obvious at the level of the price 
time series, this is not completely surprising as markets are driven by humans, who are clearly 
not time reversal invariant. In particular, the market participants remember the past, and this 
memory creates an asymmetry in the time direction. More interesting is the fact that different 
FX rates show a clear and consistent quantitative pattern. On the modelling side, the ARCH 
processes can accommodate the empirical finding, using multiscale processes with increasing 
return time horizons dt r for increasing volatility time horizons dt a . The simplest GARCH(1,1) 
process can only reproduce some asymmetry, because it contains only one time scale. On 
the other hand, all the stochastic volatility and regime switching processes are essentially time 
reversal invariant, at odds with the empirical data. This deficiency is related to the structure 
of these processes, and is therefore a key shortcoming of this class of models. Eventhough the 
various theoretical processes have very different structures, and in particular are not nested, TRI 
provides us with a strong selection criterion: a large number of processes cannot reproduce the 
stylized facts related to the time irreversibility observed in empirical time series. 
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